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INTRODUCTION 


At  the  beginning  of  this  program ,  an  evaluation  of  the  literature  in  the 
solid  propellant  extinction  field  showed  that  there  was  a  considerable  amount 
of  good  experimental  data  available  (1,2,3).  Although  some  theoretical  work 
had  also  been  done  (4),  it  was  felt  that  the  initial  efforts  would  be  most 
profitably  directed  toward  improving  the  status  of  the  theory.  Accordingly  the 
program  was  set  up  to  first,  develop  an  improved  theoretical  model;  second  to 
compare  the  mcdel  with  existing  experimental  data;  third,  to  getner  any  addi¬ 
tional  experimental  data  necessary  and  fourth,  to  expand  either  the  experimen 
tal  or  theoretical  program  as  the  results  indicated. 

THEORY 

The  theoretical  effort  was  based  on  the  premise  that  the  results  ought  to 
be  useful  to  the  design  engineer  and  to  be  useful,  the  theory  must  not  be  too 
complex  or  sophisticated.  At  least,  this  seems  to  be  true  in  the  closely  re¬ 
lated  field  of  combustion  instability  where  a  sizeable  output  of  complex 
theory  has  had  practically  no  impact  on  the  design  field.  Therefore,  a  simple 
thermal  model  was  chosen  to  represent  the  combustion  process.  This  thermal 
model  was  indicated  first  because  the  Von  Elbe  treatment  agreed  with  experi¬ 
ment  in  many  cases  and  second,  because  the  characteristic  heat  up  time  of 
the  solid  propellant  seemed  to  be  the  only  characteristic  time  of  the 
combustion  process  which  had  the  same  magnitude  as  the  extinction  time  for 
the  propellant. 

Accordingly,  the  combustion  process  was  pictures  as  follows.  The 
homogeneous,  one-dimensional  solid  propellant  with  constant  thermal  proper¬ 
ties  is  heated  to  a  critical  pressure-independent  surface  temperature  where  it 
sublimes  and  then  reacts  in  the  gaseous  phase.  The  gas  phase  processes 
are  comparatively  rapid  and  for  pressure  transients  just  fast  enough  to 


extinguish  the  combustion,  the  gas  behavior  is  dependent  only  on  the  instan¬ 
taneous  pressure.  Extinction  occurs  when  the  transient  pressure  causes  the 
propellant  surface  temperature  to  drop  below  the  critical  value. 

Appendix  1,  which  is  a  copy  of  a  paper  presented  at  the  June,  1967 
ICRPG/AIAA  Solid  Propulsion  Symposium  in  Anaheim,  California,  presents  a 
detailed  development  of  the  theory.  This  appendix  also  presents  the  theore¬ 
tical  results  in  some  detail. 

Where  possible,  the  theoretical  predictions  were  compared  with  the 
experimental  data  available  in  the  literature.  Qualitative  agreement  was 
realized  from  a  qualitative  comparison.  However,  quantitative  comparison 
could  not  be  made  because  the  data  in  the  literature  were  not  presented  in 
sufficient  detail. 

EXPERIMENT 

To  remedy  this  deficiency,  the  requisite  extinction  data  were  gathered 
experimentally.  The  combustion  chamber  was  patterned  after  that  of  Ciepluch 
(1)  and  the  method  of  initiating  the  pressure  transient  was  patterned  after  that 
used  at  the  University  of  Utah  (5).  Thus  the  propellant  would  be  ignited  in  a 
small  rocket  motor,  the  motor  would  reach  design  pressure,  and  then  a  pres¬ 
sure  transient  would  be  induced  by  rupturing  a  frangible  diaphragm  which 
sealed  an  auxiliary  nozzle. 

Both  the  experimental  procedure  and  results  are  described  in  greater 
detail  in  Appendix  1  and  Table  1  which  present  the  results  obtained  from  the 
testing  of  seven  propellant  variations.  For  the  sake  of  simplicity,  the  table 
is  based  on  the  approximation  that  the  pressure  transients  are  exponential. 

COMPARISON  OF  THEORY  AND  EXPERIMENT 

As  can  be  seen  from  Table  1,  the  theory  developed  in  this  program  agrees 
with  experiment  better  than  the  earlier  work.  This  improvement  is  accompanied 


2 


however,  by  increased  difficulty  of  application.  Also,  even  with  this  improve¬ 
ment,  the  theory  is  not  yet  adequate  for  design  pruposes.  No  comparison  was 
made  between  this  theory  and  the  more  sophisticated  theories  now  in  existence 
because  the  more  sophisticated  theories  cannot  yet  be  quantitatively  related 
to  the  experimental  data. 

ACKNOWLEDGEMENT 

Special  acknowledgement  must  be  made  to  the  personnel  cf  the  University 
of  Utah  and  Hercules  Inc.  whose  cooperation  has  eliminated  duplication  of 
effort  and  greatly  assisted  this  work. 

REFERENCES 

1.  Ciepluch,  C.C.,  "Effect  of  Rapid  Pressure  Decay  on  Solid  Propellant 
Combustion."  ARS  Journal  Vol.  31,  Page  1584,  1961. 

2.  Ciepluch,  C.C.,  "Effect  of  Composition  on  Combustion  of  Solid  Propel¬ 
lants  During  a  Rapid  Pressure  Decrease.”  NASA  Technical  Note  D-1559, 
1962. 

3.  Ciepluch,  C.C,,  "Spontaneous  Reignition  of  Previously  Extinguished 
Solid  Propellants,"  NASA  Technical  Note  D-2167,  1964. 

4.  Von  Elbe,  G.,  "Theory  of  Solid  Propellant  Ignition  and  Response  to 
Pressure  Transients."  Interagency  Solid  Propulsion  Meeting  3:96.  1963. 

5.  Donaldson,  A.B.,  "Flame  Extinction  of  Solid  Propellants,"  Unpublished 
Masters  Thesis,  University  of  Utah,  June,  1965. 


TABLE  I 


Comparison  of  Theory  and  Experiment 


JDellant* 

Initial 

pressure 

(psia) 

Experimental 
half  time** 
(sec) 

Reference  4 
Theoretical  half  time 
(sec) 

This  Theory 
half  time 
(sec) 

A-13 

165 

.0100-. 0107 

.0040 

.010-. 0107 

A-13 

410 

.006-009 

.0016 

.006-. 009 

A-13 

90 

.0105-. 012 

.0074 

.0105-. 012 

A-14 

280 

.0056-. 0085 

.00064 

.0028-. 0043 

A-14 

155 

.0057-. 0085 

.0017 

.0043-. 0058 

A-16 

320 

.0016-. 0022 

.00019 

.0008-. 0011 

A-15 

285 

.0053-. 0075 

.00093 

.001-. 0015 

A-17 

170 

.0067-. 0079 

.046 

.035-. 04 

A-18 

155 

.0062-. 0070 

.028 

.031-. 035 

125 

.0068-. 0085 

.0035 

.0023- .0041 

*  See  appendix  1  for  compositions 

**  All  of  these  runs  were  made  to  a  final  pressure  of  12.5  psia.  The  interval 
represents  the  area  separating  a  non-extinguishing  run  from  one  which 
extinguished . 
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APPENDIX  1 


DEPRESSURIZATION  INDUCED  EXTINCTION  BY 
BURNING  SOLID  PROPELLANT1 

by 

2  3  4 

M.D.  Horton  ,  P.S.  Bruno  ,  E.C.  Graesser 

ABSTRACT 

This  paper  presents  the  results  of  a  combined  theoretical  and  experimental 
study  of  the  process  in  which  solid  propellant  combustion  is  terminated  by  a 
rapid  drop  in  pressure.  Experimentally,  propellant  burning  in  a  small  rocket 
motor  was  subjected  to  a  rapid  pressure  drop  when  an  auxiliary  nozzle  was 
suddenly  opened.  Such  tests  were  run  with  varying  nozzle  sizes  for  the  aux¬ 
iliary  nozzle  and  the  boundary  between  extinction  and  non-extinction  was 
determined.  The  results  were  then  compared  to  the  predictions  of  the  theoreti¬ 
cal  model  which  was  based  on  the  assumption  that  extinction  occurs  when  the 
heat  absorption  by  the  solid  propellant  exceeds  the  heat  transfer  to  the  solid 
propellant  from  the  combustion  gas.  In  general,  the  theoretical  predictions 
agreed  well  with  both  the  experimental  results  gathered  in  this  program  and 
those  published  by  other  investigators . 


*Part  of  the  research  was  sponsored  by  AF-AFOSR  grant  No.  897-65  and  part 
by  Contract  AF  04  (694)— 1 27 ,  WS-133A,  1965  Production  Support  Program 
Task  8. 

2 

Assistant  Professor  of  Chemical  Engineering,  Brigham  Young  University. 

3 

Senior  Development  Engineer,  Hercules  Incorporated,  Magna,  Utah. 

4 

Project  Engineer  Supervisor,  Hercules  Incorporated,  Magna,  Utah. 
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INTRODUCTION 


Many  solid  propellant  applications  are  dependent  upon  the  designer's 
ability  to  cause  a  termination  of  the  propellant  combustion.  For  example,  a 
stop-start  capability  may  be  desired  for  space  missions,  and  a  termination  is 
programmed  for  the  last  stage  of  a  ballistic  missile.  Because  of  these  needs, 
engineers  have  developed  and  used  non-optimal  termination  techniques. 
Subsequent  work  in  the  area  has  been  directed  toward  understanding  extinc¬ 
tion  and  improving  this  portion  of  a  rocket's  performance. 

The  experimental  study  of  propellant  extinguishment  by  pressure  decreases 
has  been  performed  largely  in  two  types  of  apparatus.  In  one,  the  rarifaction 
tube,  a  pressure  transient  is  produced  by  bursting  a  diaphragm  at  the  end  of 
a  tube  which  contains  a  burning  strand  of  propellant.  The  resultant  pressure 
wave  then  impinges  upon  the  propellant  whose  response  to  the  transient  is 
observed.  The  other  technique  involves  firing  a  small  motor,  allowing  it  to 
reach  operating  pressure,  and  then  stepwise  increasing  the  nozzle  area  to 
cause  a  pressure  transient. 

Price ^  and  McCune^  reported  early  exploratory  tests  performed  with 

(3) 

rarifaction  tubes.  However,  Mitchell'  '  later  showed  that  such  extinction  as 
was  previously  observed  had  probably  been  caused  by  the  flow  of  cold,  inert 
gas  around  the  propellant  strand.  Donaldson ^  later  performed  an  extensive 

series  of  rarifaction  tube  tests  and  obtained  good  qualitative  results. 

(5  6  7) 

Cieplucfr  '  '  '  studied  the  extinction  of  several  propellants  by  the  use 
of  a  small  motor  having  an  auxiliary  nozzle  which  opened  after  steady  state 
combustion  was  attained.  The  reports  of  his  work  describe  the  effect  of  pro¬ 
pellant  composition  on  reignition,  and  also  the  depressurization  rate  required 
to  produce  extinction.  Reference  8  also  describes  results  of  this  nature. 
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The  theoretical  studies  performed  thus  far  have  been  only  partly  success¬ 
ful.  References  8  and  9  develop  the  transient  burning  rate  equation: 

_r  _  ,  ,  anN  dP  (1) 

r0  i+b2p2nHdt 

which  says  extinction  occurs  (r  =  0)  when 

-dP  =  b2P2n+1  (2) 

dt  Nan 

According  to  Reference  8,  N  has  a  value  of  one  or  less,  while  Reference  9 
says  N  has  a  value  of  2. 

Equation  2  provides  a  qualitative  correlation  of  the  experimental  data, 
but  a  better  theoretical  description  is  required  for  two  reasons.  One  reason 
is  that,  as  Figure  1  show,  the  correlation  is  not  adequate  for  general  engi¬ 
neering  use.  Another  reason  is  that  Equation  2  does  not  consider  the  effect 

(7) 

of  total  pressure  drop  which  experiment  has  shown  to  be  important. 

The  object  of  this  paper  is  to  present  the  results  of  a  study  in  which  an 
improved  theoretical  model  explaining  extinction  was  developed.  Also  pre¬ 
sented  are  experimental  extinction  data  for  propellant  compositions  not  pre¬ 
viously  tested. 


EXTINCTION  MODEI, 

To  derive  a  rigorous  theory  describing  extinction,  one  would  develop  the 
time  dependent  heat  conduction,  diffusion,  and  kinetic  equations  for  both  the 
gas  and  solid  phases  of  the  combustion  zone,  and  then  solve  them  simulta¬ 
neously.  Besides  being  extremely  difficult,  this  task  would,  upon  completion, 
yield  a  very  complex  solution.  In  fact,  the  solution  would  be  so  complex  and 
contain  so  many  parameters  whose  values  were  unknown  that  such  a  solution 
would  have  little  practical  value. 
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The  task  could,  perhaps,  be  simplified  by  the  use  of  perturbational 
analysis,  except  that  here  again  difficulties  intervene.  During  the  pressure 
transient  being  considered,  many  pertinent  variables  change  by  a  factor  of 
ten  or  greater.  This  being  so,  one  must  be  concerned  not  with  a  perturbation 
of  the  initial  conditions,  but  rather  with  the  entire  transient  of  interest. 

Fortunately,  the  problem  can  be  somewhat  simplified  if  consideration  is 
given  to  the  various  characteristic  times  for  the  combustion  process.  Reference 
10  has  shown  that  the  characteristic  times  for  the  gaseous  portion  of  the  com¬ 
bustion  zone  are  considerably  smaller  than  the  characteristic  time  associated 
with  heat  conduction  in  the  solid  propellant.  Accordingly,  the  heat-up  of  the 
solid  propellant  is  the  slowest  in  the  series  of  events  that  transforms  the 
solid  propellant  into  a  hot  combustion  gas.  Because  the  heating  of  the  solid 
is  the  rate  limiting  (i.e.  slowest)  step,  the  minimum  pressure  transient  which 
will  terminate  combustion  is  the  one  which  is  just  fast  enough  to  arrest  this 
rate-limiting  step. 

In  order  to  derive  a  thermal  extinction  criterion,  it  is  first  convenient  to 
describe  the  steady  state  thermal  profile  in  the  solid  propellant.  The  heat 
conduction  equation  in  the  solid  propellant  is: 

,  32T  _  r  3  T  _  3  T  (3) 

3  2  ri  3  x  3  t 

3  X 

if  it  is  assumed  that: 

(1)  the  heat  conduction  is  one-dimensional 

(2)  there  are  no  reactions  beneath  the  propellant  surface 

(3)  the  thermal  properties  of  the  propellant  are  constant 

(4)  the  propellant  is  homogenous  and  semi-infinite,  and 

(5)  the  coordinate  system  is  fixed  with  respect  to  the  propellant 
surface 
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At  steady  state  3  T  =  0  and  the  temperature  distribution  in  the  solid  propel- 
3  t 

pant  is  (see  Figure  2): 


T  -  T 


s  -  T 


o _  ,rxN 

-  exp  (— ) 


(4) 


Now,  it  is  further  assumed  that  the  propellant  sublimes  with  a  pressure- 
independent  heat  of  sublimation  when  it  is  heated  to  a  certain  critical, 
pressure-independent  temperature  (Ts).  Certainly,  this  is  not  an  exact  des¬ 
cription  of  the  surface  decomposition,  but  it  may  be  an  adequate  representa¬ 
tion  as  will  be  determined  by  a  comparison  of  the  theoretical  predictions  and 
experimental  data . 

The  consequence  of  these  assumptions  is  that  any  time  the  surface 
temperature  drops  below  T£,  extinction  has  occurred.  Furthermore,  the  only 
way  the  surface  temperature  can  drop  is  the  for  the  heat  conduction  into  the 
propellant  to  be  greater  than  the  heat  conduction  from  the  gaseous  flame  to 
the  propellant  surface. 

3  T  l 

The  heat  conduction  into  the  solid  is  k  —  _  n  whereas  it  is  more 

difficult  to  describe  the  conduction  to  the  surface.  To  facilitate  this  descrip¬ 
tion,  another  approximation  is  made.  This  approximation  is  that  the  rate  of 
heat  conduction  from  the  gas  flame  to  the  solid  surface  is  characterized  only 
by  the  pressure.  That  is,  at  a  given  pressure,  the  heat  transfer  is  the  same 
during  a  pressure  transient  as  it  would  be  if  the  pressure  were  constant.  To 
some  degree,  the  approximation  is  justified  because  the  characteristic  times 
for  the  gas  are  small  compared  with  the  characteristic  time  of  the  solid.  ^ 

This  means  that  the  gas  flame  responds  rapidly  enough  that  its  behavior  is  not 
time-dependent,  but  only  pressure  dependent.  Again,  the  ultimate  justification 
must  come  from  satisfactory  agreement  between  experiment  and  theory.  The 
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heat  conducted  into  the  solid  from  the  surface  is  equal  to  that  conducted  to 

the  surface  minus  the  heat  required  to  vaporize  the  solid.  With  the  above 

3  T I 

assumptions  and  approximations,  these  respective  terms  are  k- —  _  n 

d  X  |  X*  —  \j  t 

rpi  +  r  pc  (T  -  T  ),  and  r.  pA  .  The  equation  which  relates  these  terms  is: 
so  i 


fr|x  =  0,  =p>(r-ri)+r|>cCrs'To) 


(5) 


At  this  point  Equation  3  describes  the  heat  transfer  in  the  solid. 

Equation  4  is  a  boundary  condition  describing  the  initial  temperature  profile 
when  the  pressure  transient  is  induced,  and  Equation  5  describes  the  transient 
heat  transfer  to  the  solid  surface.  Additional  boundary  conditions  are  that  T 

b 

is  a  constant  and  the  temperature  at  minus  infinity  is  T  . 

Eor  the  sake  of  generality,  these  equations  were  non-dimensionalized  by 
the  use  of  the  following  definitions,  many  of  which  were  used  in  Reference  8: 

T  -  T 


0  = 


T  -  T 
s  o 


(6) 


a 

x  - 


R  =  - 
ro 


r, 

'i  =  r. 


R,  =  r1- 


F  = 


c<Ts-V 


The  substitution  of  these  dimensionless  variables  converted  the  heat- 
transfer  equations  and  boundary  conditions  to: 
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3  2  39  _  39 

ay2  ‘  Ki  3Y  3t 

-fH  Y-  0  -  tl  +  F)  R  -  RjF  (7) 


9  -  0  at  Y  -  -® 


9  -  exp  Y  at  t  =■  0 


and 


9  =  1  at  Y  =•  0 


Provided  that  the  relationship  between  R  and  x(p  and  t)  were  known,  one 
could,  in  principle,  find  an  exact  solution  to  this  set  of  equations.  However, 
even  with  a  simple  specification  of  R  versus  x  ,  the  authors  could  not  find  an 
exact  solution  and  resorted  to  numerical  techniques.  By  the  application  of  the 
backward  method  of  differences,  a  tridiagonal  system  of  difference  equations 
was  developed  and  then  numerically  solved  with  an  IBM  7040  computer  for  two 
types  of  Rvs.x  relationships. 

The  first  type  was  chosen  because  it  was  a  fair  approximation  of  the  rela¬ 
tionship  encountered  experimentally.  This  form  was  based  on  the  assumption 
that  the  pressure  decay  in  the  combustion  system  was  described  by  the 
exponential  equation: 


P  -  P 
_ a_ 

P  -  P 
o  a 


=  exp  ( 


-0.693t 


(8) 


while  the  steady  state  burning  rate  of  the  propellant  could  be  described  as: 


r 


bP 


n 


(9) 


Then,  the  equation 


(— )n 
(po 


(10) 


11 


was  used  along  with  Equations  8  and  10  to  obtain 

«-  ca-  exp  (~-:p~  >  +  /] 

o  %  ro 

Equations  7  and  11  were  then  solved  numerically  for  the  critical  values 
of  which  represented  the  boundary  condition  between  fast  decays 
(small  which  theoretically  caused  extinction  (R^  =0)  and  slower  decays 
which  did  not.  Figures  3  and  4  summarize  the  results. 

In  many  cases  of  interest,  either  or  both  of  Equations  8  and  9  did  not 
adequately  describe  the  transient.  When  this  was  so.  Equation  11  was  not 
used  as  a  boundary  condition  but  instead,  the  computer  program  used  a  data 
table  of  R  versus  t  as  prepared  from  the  strand  burning  rate  curve  and  the 
pressure  decay  record  of  a  particular  extinction  test.  The  numerical  solution 
would  then  determine  whether  or  not  the  propellant  should  have  extinguished 
theoretically.  The  results  of  such  tests  were  examined  by  this  procedure  and 
are  presented  later  in  the  paper. 

EXTINCTION  TESTS 

Although  much  experimental  extinction  data  were  available  in  the  litera¬ 
ture,  they  did  not  meet  the  needs  of  this  program  for  two  reasons.  First, 
because  propellant  compositions  not  previously  tested  were  of  interest  and 
second,  because  the  data  as  presented  did  not  allow  a  comparison  between 
experiment  and  theory.  That  is,  the  pressure  transients  in  the  literature 
were  non-exponential,  so  Equation  11  was  not  a  valid  boundary  condition. 
Also,  the  experimental  decays  were  rarely  presented  in  a  sufficient  detail 
to  permit  the  determination  and  use  of  a  tabular  boundary  condition.  One 
experimenter  was  generous  enough  to  furnish  the  authors  with  his  detailed 
experimental  results.  However,  the  data  were  originally  taken  for  dif¬ 
ferent  purposes  and  were  not  detailed  enough  to  allow  a  direct  comparison 
with  theory. 
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Therefore,  the  propellant  compositions  listed  in  Table  I  were  tested  to 
determine  the  conditions  required  to  extinguish  the  propellant  combustion. 
This  testing  was  done  in  motors  which  were  ignited,  allowed  to  come  to 
operating  pressure,  and  then  subjected  to  a  pressure  transient  by  instantly 
opening  an  auxiliary  nozzle.  The  use  of  several  different  auxiliary  nozzle 
sizes  established  the  boundary  between  the  transient  which  extinguished 
combustion  and  the  transient  which  did  not.  Studied  in  the  program  were  the 
effect  of  propellant  location,  propellant  composition,  initial  pressure,  and 
final  pressure.  Table  2  presents  a  summary  of  the  test  results. 

As  Reference  6  has  shown,  small  changes  in  composition  can  cause  large 
changes  in  the  extinguishment  characteristics  of  a  propellant.  Other  vari¬ 
ables  had  less  effect  but  were  still  significant. 

DETERMINATION  OF  F 

The  reader  may  have  noticed  that  all  but  one  of  the  parameters  used  in 

the  theory  can  be  readily  determined  from  routine  laboratory  tests .  Recent 
(12) 

work  has  also  made  it  possible  to  determine  the  value  of  F.  Briefly,  the 
value  of  F  is  determined  as  follows.  A  rarifacticn  tube  is  used  to  determine 
the  total  pressure  drop  required  to  extinguish  a  burning  propellant  strand  as 
a  function  of  depressurization  rate.  This  curve  is  then  extrapolated  beyond 
the  experimental  limit  to  fine  the  required  pressure  drop  (PQ  to  P^)  for  an  in¬ 
stantaneous  drop  in  pressure.  This  parameter  is  then  related  to  F  oy  noting 
that,  according  to  the  assumptions  made  earlier,  an  instantaneous  pressure 
drop  will  cause  the  heat  flux  from  the  gas  phase  to  instantaneously  readjust 
while  the  thermal  profile  in  the  solid  is  unchanged.  The  extinction  condition 
is  represented  by  the  point  where  the  heat  conduction  to  the  surface 

rpA  +  rpc  (T  -  T  )  is  just  equal  to  the  conduction  away  from  the  surface 
s  c 

ro  pc  (Ts  -  V  and 
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(12) 


rp[x+c(Ix-T0)]=r0Pc  tts  -  T„) 


Equation  9  is  then  used  as 


r  =  bP 


n 


(13) 


and 


r  =  bP 
o  o 


n 


then  combined  with  Equation  12  to  give: 


{P°}  £  c  ( T  -  T  )  +11 

X  so 


(14) 


which,  after  the  substitution  of  F  and  rearrangement,  becomes 


P  n 

F  =  (~)  -1 
*f 


(15) 


For  several  of  the  propellants  tested  in  the  extinction  program,  the  values 
of  F  and  the  parameters  used  in  their  calculation  are  shown  in  Table  III. 


COMPARISON  OF  THEORY  AND  EXPERIMENT 
Even  though  the  theoretical  model  could  not  quantitatively  be  compared 
with  the  experimental  data  available  in  the  literature,  some  qualitative  com¬ 
parisons  were  made.  Ciepluch  (Figure  3  in  Reference  6)  shows  a  graph  of 
the  maximum  "half-decay  time"  (t^  ),  which  will  cause  extinction  as  a 

function  of  initial  pressure.  From  the  definition  of  t  ,  we  find  that 

2 

T  _  ro  t  (16) 

hah 

Now,  to  compare  theory  and  experiment,  the  rather  poor  assumption  was  made 
that  the  experimental  pressure  decay  is  exponential.  Then,  for  reasonable 
values  of  F  (see  Tables  2  and  4),  Figures  3  and  4  show  thatx^  is  a  constant 
over  the  pressure  range  being  considered  (n  =  .34,  P  /P  =  .001).  As  figure 

3  O 

S  shows,  the  experimental  data  have  about  the  same  values  as  the  theory 
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(F  =  0.1),  but  the  experimental  data  show  a  slight  slope.  Accordingly,  the 
model  and  experimental  data  show  qualitative  agreement  which  is  about  all 
that  could  be  expected  in  this  case. 

Reference  6  presents  data  for  the  critical  half  decay  time  as  a  function 
of  ambient  pressure  for  an  aluminized  propellant.  Figure  6  shows  a  compari¬ 
son  between  the  theory  and  the  experimental  results  and  again  qualitative 
agreement  is  observed. 

In  one  case,  an  illustration  in  the  literature'^  presented  exact  pressure¬ 
time  history  during  an  extinction  test.  This  history  was  used  to  prepare  a 
table  of  R  vs.  t  for  computer  input.  Then,  with  an  assumed  representative 
value  of  F  =  0.10,  the  theory  correctly  predicted  extinction.  Because  an 
exact  value  of  F  was  unknown,  this  result  represents  qualitative  agreement 
between  experiment  and  theory. 

Similarly,  the  F  values  for  several  of  the  propellants  tested  in  this  pro¬ 
gram  were  unknown.  However,  the  use  of  an  assumed  value  (0.14)  and 
tabular  input  of  the  boundary  condition  yielded  excellent  agreement  between 
the  experimental  results  and  the  theoretical  predictions  with  agreement  being 
obtained  in  fourteen  out  of  fifteen  comparisons.  This  again  constitutes  good 
qualitative  agreement. 

For  those  propellants  whose  F  values  were  known  (See  Table  III) ,  use  of 
tabular  input  of  the  boundary  condition  provided  a  severe  test  for  the  theoreti¬ 
cal  model.  These  seven  propellant  compositions  were  tested  and  the  theory 
correctly  predicted  extinction  for  18  out  of  25  tests.  To  determine  how  serious 
the  error  in  the  model  was,  the  seven  tests  showing  disagreement  were  further 
examined.  Faster  or  slower  decays  were  synthetically  generated  and  used  to 
determine  the  extinction-non-extinction  boundary.  By  this  process  it  was 
found  that  to  predict  extinction  the  model  required  depressurization  rates 
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which  were  in  error  by  as  much  as  a  factor  of  five.  For  example,  to  predict 
extinction  for  run  119,  the  model  needed  a  depressurization  rate  at  least 
three  times  as  fast  as  the  experimentally  observed  value. 

It  is  appropriate  that  this  theory  be  compared  with  that  presented  in 
References  8  and  9.  The  perturbational  approach  used  in  those  references 
yields  a  criterion  which  is  very  simple  to  use  and  which  disagreed  by  no 
more  than  a  factor  of  ten  with  the  experimental  results  determined  in  this  pro¬ 
gram.  However,  that  approach  wrongly  predicts  no  relation  between  extinction 
and  ambient  pressure.  The  model  presented  in  this  paper  requires  complex 
numerical  calculation  but  seems  to  be  correct  within  a  factor  of  five  and  cor¬ 
rectly  predicts  trends.  It  would  seem  then,  that  the  older  theory  is  not  as 
good  but  is  easier  to  use. 

One  factor  which  was  examined  experimentally  was  the  effect  of  propellant 
orientation  in  the  combustion  chamber.  It  was  found  that  an  end  burner  was 
slightly  easier  (25  percent  smaller  depressurization  rate)  to  extinguish  than 
a  single  slab  subjected  to  parallel  flow  which  in  turn  was  slightly  easier 
(20  percent)  to  extinguish  than  opposed  slabs  subjected  to  parallel  flow. 

The  practical  implications  of  Figures  3  and  4  are  many.  According  to  the 

theory,  anything  that  lowers  the  effective  propellant  surface  temperature  or 

raises  its  "heat  of  sublimation"  will  make  the  propellant  easier  to  extinguish. 

For  low  values  of  P  /  P  ,  a  small  increase  in  either  chamber  pressure  or 
9  O 

pressure  exponent  can  make  a  motor  much  easier  to  extinguish.  Also,  under 
any  conditions,  very  low  exponent  propellants  are  difficult  to  extinguish. 
Unfortunately,  the  experimental  data  required  to  test  these  predictions  are 
not  available. 

As  the  model  was  used,  it  was  necessary  to  know  the  experimental  pres¬ 
sure  decay  to  see  if  the  theory  would  predict  extinction.  The  theory  would  be 
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much  more  useful  if  it  were  coupled  to  a  mass  balance  on  the  combustion 
chamber  so  that  the  computer  would  generate  the  pressure  transient  equation. 
While  such  a  coupling  was  accomplished  in  this  program,  an  overly  simple 
mass  balance  was  used  and  high  quality  results  were  not  obtained. 

DISCUSSION 

At  this  point,  consider  how  good  the  theory  could  be  expected  to  be  and 
then  how  good  it  actually  is.  The  theoretical  model  is  relatively  simple  and 
unsophisticated.  The  solid  is  not  homogeneous,  the  thermal  properties  are 
not  constant,  the  surface  temperature  is  not  constant,  subsurface  reactions 
almost  surely  occur,  the  use  of  the  term  X  represents  a  drastic  simplification 
of  the  surface  decomposition,  and  the  representation  of  the  gas  phase  heat 
transfer  may  be  a  poor  approximation.  V/ith  these  objections  in  mind,  one 
would  not  expect  more  than  qualitative  agreement  between  theory  and  experi¬ 
ment. 

Now,  the  comparisons  made  showed  qualitative  agreement  in  all  cases 
and  quantitative  agreement  in  most  cases.  This  is  better  than  one  would 
expect,  considering  the  simplicity  of  the  theory  and  the  wide  range  or  pro¬ 
pellant  compositions  considered.  The  theory  then,  consists  of  a  fairly 
accurate  thermal  representation  of  the  extinction  process.  Because  all  of 
the  parameters  used  can  be  determined  experimentally,  no  "guess"  factors 
are  left  as  uncertainties.  Accordingly,  the  theory  and  experimental  method 
ought  to  be  useful  to  anyone  concerned  with  practical  extinction  problems. 
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NOMENCLATURE 

Exponential  pressure  decay  exponent 

Constant  in  Vielle  steady  state  burning  rate  equation 

Propellant  heat  capacity 

Ratio  of  propellant  heat  of  vaporization  to  heat  required  to 
raise  propellant  to  surface  temperature 

Propellant  thermal  conductivity 

Pressure  exponent  in  Vielle  burning  rate  equation 

Coefficient  in  Von  Elbe  and  Aerojet  Extinction  Theory 

Instantaneous  chamber  pressure 

Ambient  pressure 

Final  pressure  in  rarifaction  tube 

Initial  steady  state  chamber  pressure 

Steady  state  burning  rate  corresponding  to  P 
Transient  burning  rate 

Steady  state  burning  rate  corresponding  to  Po 

Dimensionless  steady  state  burning  rate  (  — ) 

r.  0 

Dimensionless  transient  burning  rate  (~) 

ro 

Time 

Largest  time  in  which  charnV-  Pressure  could  decrease  by  50% 
and  still  extinguish  propr 1 

Temperature  of  propellant 

Initial  propellant  temperature 

Surface  temperature  of  burning  propellant 

Distance  from  propellant  surface 

18 


Dimensionless  distance  (rQx/ a  ) 

Propellant  thermal  diffusivity 
Propellant  heat  of  sublimation 
Propellant  density 

(L= 

Dimensionless  temperature  Tg 

2 

Dimensionless  time  (r  t/ct  ) 
Dimensionless  half  decay  time  ^ 
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PROPELLANT  COMPOSITIONS 


XJ 

XJ 

< 


© 

3 

Cm 


<0 


© 

Si 


XJ 

•H 


«  c 
C  o 
(0  t2 
to 

gl 

-a 


a 

« 

\ 

Q. 

< 


XJ 

0) 


XJ 

O 

£ 


to 

to 

o 


K>  © 
C  to 
O  io 
— i  o 

C  ' 
<D  ® 

>s 

G  -3 
°  ° 
O  XJ 


A 


5 


a, 

< 


XJ 

© 


>.JS 


s-q 

2  o 

CO  XJ 


CQ 


< 

< 

< 

55 

55 

55 

O 

CO 

o 

CM 

CM 

A 

2 

a 

* 

< 


XJ 

© 


XJ 

O 

£ 


xj 

o 

£ 

>  © 
4-J  CO 
CO  (0 
(0  XJ 
O  I 


co 

© 

O 


o 


IZ 


a  a) 
c  co 
o  to 

— «  o 

s© 

cl 

Oxj 


* 

* 


© 


© 


£ 

o 

u 

XI 

O 


£ 

o 

u 


© 

a 

a 

O  X3 

o  © 
$$  xj 

r-H  10 


x: 

X) 

^ 3 

o 

© 

© 

t. 

X3 

XI 

© 

XI 

XJ 

a 

© 

© 

a 

o 

XJ 

Cm 

m-4 

Cm 

**H 

o 

© 

X) 

>-J 

55 

XJ 

XJ 

55 

55 

fH 

© 

f-H 

H 

2 

5 

2 

53 

£ 

< 

< 

< 

< 

< 

< 

CQ 

CQ 

CQ 

CQ 

CQ 

CQ 

a, 

a. 

CM 

a. 

CL. 

a, 

55 

55 

55 

55 

55 

55 

•M* 

•sr 

M* 

■M* 

M* 

CM 

CM 

CM 

CM 

CM 

CM 

© 

© 

© 

© 

© 

© 

•H 

•H 

■H 

■H 

■H 

co 

co 

co 

CO 

CO 

co 

o 

o 

O 

o 

o 

o 

a 

a 

a 

a 

a 

a 

£ 

£ 

s 

E 

£ 

£ 

o 

o 

o 

o 

0 

o 

O 

O 

O 

O 

O 

O 

co 

•M" 

CO 

CO 

00 

•H 

| 

| 

. — i 

1 

»H 

t 

7* 

< 

< 

< 

< 

< 

< 

a. 

a 

CL, 

a, 

a. 

cu 

a 

< 

< 

< 

< 

< 

< 

< 

a 

a 

a 

a 

a 

3L 

XJ 

o 

CO 

o 

to 

o 

LO 

co  © 

00 

f-H 

00 

«-H 

00 

fH 

©  Z 

55 

55 

55 

55 

55 

55 

°  © 
00  u 

CO 

co 

CO 

co 

co 

co 

r*. 

tN. 

c^. 

* 

* 

* 


CQ 

a, 


55 

o 

CM 


© 


co 

O 

a 

£ 

o 

O 


O 


© 


w 

--  XJ 


c 

o 


o 

© 


o 


2  & 
© 
u 

o 


g  3 
©  © 


JC 

o 
1-1 

©  _ 

a  V 


■  i 
©  © 


c  c 
©  © 


£ 

3 


_  XJ 
2  © 


c 

o 

£ 

£ 

© 


3 

XI 

>. 


3 

XJ 


o 

a 


o 

a 


2  «2 


co 


a, 

< 

* 


< 

CQ  CQ 


a. 

* 

* 


O. 

* 

* 

* 


!  < 


in 

to 

t£> 

to 

to 

to 

to 

to 

to 

to 

to 

to 

to 

to 

to 

co 

CO 

CO 

CO 

co 

eg 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

co 

CO 

CO 

CO 

CM 

CM 

CM 

CM 

CM 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

XT 

"M4 

<3* 

TT 

TT 

00 

00 

CO 

CO 

co 

rH 

H 

i—i 

rH 

r—i 

i—i 

rH 

r— H 

rH 

r—i 

rH 

rH 

rH 

rH 

r-H 

rH 

r— < 

rH 

i—i 

»— H 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

(fl 

W 

(/) 

(/) 

0) 

CO 

CO 

CO 

CO 

(0 

Q) 

<s 

0) 

0 

a) 

a) 

0) 

o 

Q 

<D 

0 

0 

0) 

0 

o 

o 

o 

<s 

o 

o 

>1 

>i 

is 

>, 

>1 

>i 

>1 

c 

C 

>1 

>1 

c 

>1 

c 

c 

c 

c 

>1 

c 

c 

CO 

CO 

CO 

c 

CO  CO 

CO 

(0 

CO 

CO 

CO 

0) 

0 

0 

o 

rlj 

0  0 

0 

o 

o 

0 

0 

o 

0 

o 

o 

O 

o 

© 

o 

o 

>1 

>1 

>1 

c 

u 

a) 

>1  >1 

>1 

c 

G 

>1 

>1 

c 

>1 

c 

G 

c 

G 

c 

c 

o 

jf«0 

U  03 


oo 

CO 

o 

o 

CO 

o 

r-H 

co 

CO 

CO 

V 

T 

10 

o 

co 

00 

00 

CM 

■"S' 

CO 

o 

o 

o 

03 

O 

o 

o 

o 

o 

CM 

CM 

v 

rH 

H 

r—i 

CM 

co 

CO 

r—i 

rH 

r—i 

r—i 

O 

O 

• 

o 

• 

o 

o 

• 

o 

o 

• 

O 

• 

o 

o 

« 

o 

• 

o 

o 

O 

• 

o 

• 

o 

• 

o 

• 

o 

• 

o 

• 

o 

• 

O 

• 

—4 

CO 

•*r 

rH 

•M* 

00 

00 

CO 

in 

00 

o 

co 

to 

CM 

C' 

to 

CO 

CO 

03 

CM 

to 

m 

03 

o 

rH 

o 

■M* 

M* 

Tl* 

CM 

CM 

co 

• 

-M’ 

• 

CM 

• 

■M* 

CO 

• 

CM 

• 

co 

• 

CM 

• 

CO 

• 

in 

• 

• 

CO 

• 

•H 

• 

r—i 

• 

rH 

• 

r—i 

• 

CM 

• 

rH 

m 

in 

m 

m 

o 

o 

m 

to 

in 

in 

m 

m 

in 

m 

tn 

• 

• 

• 

• 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

• 

• 

• 

• 

• 

• 

o 

CM 

CM 

CM 

• 

• 

• 

• 

♦ 

• 

• 

• 

• 

CM 

CM 

CM 

CM 

CM 

CM 

rH 

rH 

•H 

o 

o 

c* 

o 

o 

o 

o 

o 

o 

o 

r—i 

r—i 

r—i 

r—i 

r—i 

rH 

o 

00 

CO 

to 

to 

a> 

(T) 

o 

o 

CM 

o 

t>~ 

m 

VP 

C- 

CO 

CO 

o 

X 

03 

co 

CO 

m 

o 

o 

o 

V 

CM 

r—i 

03 

CO 

CO 

co 

00 

to 

CO 

m 

m 

rH 

-M* 

CO 

rH 

CO 

CM 

CO 

•M* 

CO 

CM 

rH 

r—i 

r—i 

co 

*  *  *  * 
*  *  *  * 
CO  03  03  O 

CM  CO 


*  * 
*  * 
co  in 
cm  co 


*  *  *  *  * 

*  *  *  *  * 

N  l/}  CO  ^ 

co  to  in  co  co 


o  o  o  o 


cm  cm  m  h  HHNn22:110U,<91®**** 

■^OOOO  0000E"'E‘,00000Ol0u:>00 

1 — %  Ht  > — i  I— I  Hi  H- 1  Hi  > — 1  Cm  Cm  H-i  t — 1  H— 1  I  I  > — 1  •— I  r—t  H  M 

IS 

ij  CO  CO  CO  CO 

g)  iM  iH  iM  iH 

^  <1  <1  <1  <;  <!4!<!<<!<!cqOOOO<|<|i4|4J 

£ 

o 


EXPERIMENTAL  RESULTS 


CO 

co 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

co 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

C". 

•  * 

o 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

<-M 

CM 

CM 

3  ; 

fl) 

o 

o 

o 

O 

Q 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

; 

o 

o 

o 

o 

O 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

J  ‘ 

\ 

c 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

CM  . 

C 

o 

o 

o 

o 

O 

o 

CD 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

; 

i — I  «-H  H  i—t  i— t  r— I  *— I  < — t  LO  LO  <H  Hi-iHfHHHiHHH 

ooooooooooooooooooooo 


M 

W 

OT 

ui 

W 

CO 

U) 

Ul 

Ul 

Ul 

Ul 

<D 

o 

0) 

o 

eu 

o 

d) 

0) 

O 

O 

o 

o 

a) 

CD 

CD 

0) 

CD 

a> 

o 

0 

Is 

c 

is 

c 

c 

is 

is 

C 

C 

c 

c 

is 

Is 

is 

>* 

is 

is 

c 

c 

a 

JC  ® 
u  ui 


co 

co 

co 

CO 

co 

CO 

co 

co 

co 

CO 

CO 

l  ? 

CD 

O 

d> 

O 

CD 

o 

a) 

o 

a) 

o 

CD 

O 

o 

Q 

a> 

o 

a> 

0) 

O 

CD 

O 

is 

c 

is 

c 

is 

c 

is 

c 

is 

c 

is 

c 

c 

is 

is 

c 

is 

c 

>. 

C 

t—4 

to 

CO 

Cs 

to 

CO 

to 

CO 

CM 

CO 

CO 

ai 

CS 

CM 

o 

•H 

CO 

in 

co 

CO 

00 

to 

00 

t-4 

CM 

in 

ts 

t~~ 

to 

CO 

ts 

to 

ts 

o 

CO 

00 

o 

S 

o 

o 

o 

o 

o 

o 

o 

c 

o 

o 

o 

o 

o 

o 

t—4 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

ffl 

to 

o\ 

CO 

co 

CD 

o 

o 

00 

00 

CO 

t—4  ?*s> 

in 

co 

m 

ID 

o 

CM 

r—4 

o 

CM 

CM 

CM 

>tr 

Cs 

00 

00 

in 

r-v 

m 

ID 

CM 

00 

o 

cv 

co  ,1 

00 

H  • 

CM 

•H 

CM 

CM 

M*  CO 

CD 

ID 

co 

co 

t—4 

r-H 

t—4 

t—4 

•—4 

t—4 

c 

• 

• 

• 

• 

• 

• 

•  • 

• 

• 

• 

• 

+ 

• 

• 

• 

• 

• 

« 

• 

• 

m 

to 

m 

in 

in 

m 

in 

in 

in 

u> 

in 

in 

in 

in 

to 

m 

in 

m 

m 

in 

l/> 

j 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

t—4 

t—4 

r—4 

t—4 

pH 

pH 

r-H 

pH 

pH 

pH 

t—4 

r—4 

t—4 

pH 

pH 

r—4 

pH 

t—4 

r  i 

p  i 

si 

o 

CM 

CO 

CD 

00 

00 

in 

o 

CD 

o 

to 

ID 

TT 

o 

in 

o 

00 

o 

to 

"J 

•M" 

OT 

CD 

00 

m 

M* 

o 

ID 

r—4 

CM 

ts 

at 

CD 

N 

■tr 

CM 

CD 

CM 

CO 

ri 

■tr 

CM 

fH 

r—4 

CO 

CM 

00 

CO 

CM 

CM 

t—4 

pH 

pH 

pH 

CM 

CM 

CM 

t—4 

pH 

y 

*  -It 

*.fc****************.leO>a> 

oJHNT?'a,iflO>0^,<OWlOO,-<tONOOO'-t^N 

NNC'lNCOnW^'^'ViflifllDIOWlOIOSNiHrH 


7  <T  ^  V  <0  (O  U) 


N  N  CO  CO  R 


I  I  ■  !  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  ,R 


*end  burner 
**opposed  slabs 


TABLE  III 


RARIFACTION  TUBE  RESULTS 


Propellant 

Vpf 

n 

F 

A- 13 

1.265 

0.51 

0.13 

A- 14 

1.29 

0.61 

0.17 

- 

A- 15 

1.87 

0.53 

0.40 

A- 16 

2.24 

0.56 

0.57 

A- 17 

1.29 

0.42 

0.11 

A- 18 

1.204 

0.57 

0.11 

G 

1.306 

0.58 

0.17 

Figure  1  COMPARISON  BETWEEN  THEORY  AND 
EXPERIMENTAL  EXTINCTION  CRITERION 
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TEMPERATURE  PROFILE  IN  SOLID  PROPELLANT 


Figure  3  THE  MAXIMUM  DIMENSIONLESS  TIME  IN  WHICH 
THE  CHAMBER  PRESSURE  CAN  EXPONENTIALLY  DECAY  TO  HALF 
THE  VALUE  AND  STILL  THEORETICALLY  CAUSE  THE  PROPELLANT 
COMBUSTION  TO  EXTINGUISH  FOR  AN  AMBIENT  PRESSURE  OF 


400  800  1200 


PQ  (  Initial  chamber  pressure,  psia  ) 


Figure  5  THE  MAXIMUM  DIMENSIONLESS  TIME  IN  WHICH 
THE  CHAMBER  PRESSURE  CAN  DECAY  TO  HALF  THE  INITIAL 
VALUE  AND  STILL  EXTINGUISH  THE  PROPELLANT  COMBUSTION. 
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